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ABSTRACT 

Main aim of this paper is the first detailed analysis of multiple system V2083 Cyg 
and to reveal its basic physical properties. The system was studied by method of the 
light and radial velocity curves analysis, together with the interferometric data of the 
visual pair obtained during a last century. There was found that the close subsystem 
contains two very similar stars of spectral type A7-8. Moreover, the third body is 
orbiting around this pair with period of about 177 years. Due to the discrepancy of 
total mass as derived from two methods, there arises that the third body is maybe also 
a binary, or some object with lower luminosity but higher mass than normal main- 
sequence star. Another explanation is that the Hipparcos value of parallax is incorrect 
and the system is much closer to the Sun. 

Key words: binaries: eclipsing - binaries: visual - stars: fundamental parameters - 
stars: individual: V2083 Cyg. 



1 INTRODUCTION 

The eclipsing binaries as members of more complex multiple 
systems can provide us important information about their 
physical properties, as derived from different methods. This 
is the case of V2083 Cyg, which is the system, where the 
close components form an eclipsing binary, and the third 
distant body orbiting the close pair is detected as a visual 
component. Thanks to the combined analysis we are able to 
derive the radii, masses and evolutionary status of the close 
components and also some properties of the distant one. 
Such systems are still very rare and mostly lie relatively close 
to the Solar system. Nowadays, there are known only 33 
such systems, where a close eclipsing binary is a member of 
a wide visual binary and we know both orbits, their mutual 
inclinations, ratio of periods, etc. Such a unique systems are 
the most suitable ones for studies of dynamical effects, the 
short and long-t erm evolution of the orbits, etc. (see e.g. 
ISoderhielmlll975l '). 



2 THE SYSTEM V2083 CYG 

The system V2083 Cyg (= HD 184242 = HIP 96011, RA 
19'' 31™ 16.36" DE +47° 28' 52.24", Vma=o = 6.86 mag) is 
an Algol-type eclipsing binary with its orbital period about 
1.87 day. It is also a primary component of a visual double 



star designated as WDS J19 313-I-4729 in the W ashington 
Double Star Catalog (WD^. iMason et al.llJOoil '). The sec- 
ondary component of this double star is about 220 mas dis- 
tant and is a bit fainter. On the other hand, the magnitude 
difference is not very certain, because different authors list 
different values. The WDS catalogue itself gives 7.50 and 
7.93 mag for both components. 

The system is rather neglected one and there are only 
a few papers published. It was discovered as an eclip sing bi- 
nary from the Hipparcos data (jPerrvman et al.lll997^ ■ which 
also reveal that the light curve (hereafter LC) shows two 
similar minima and a classical feature of an Algol-type star. 

Spectral typ e of t he s ystem is not known very pre- 
cisely nowadays. lAbJ (|l985i ) presented spectral classifica- 
tion of the whole AB system as Am (K/H/M=A3/A8/A9), 
iRenson et al.l (|l991h gives a composite spectral type as 
A3-A9, while the s pectra l type A3 was p resented by 
ICannon fc PickerinS (|l918f ) , lOchsenbeinI (|l98(j) , and many 
others. This could indicate that the combined spectrum 
is composed from components of slightly different spectral 
types. The photometry of V2083 Cyg obtained from the Hip- 
parcos mission giv es a color inde x B — V = 0.279 mag (indi- 
cating sp.type A9. |Poppeilll980l ). while the infrared J — H 
and H — K indices, which are less influenced by interstellar 
reddening, as derived fr om the 2MA SS survey give spectral 
types about A4 and A7. ICoxl l|2000l '). 
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Figure 1. Light curves in B (blue), V (green), and R (red) filters 
for V2083 Cyg, the solid fine represents the final fit (see the text). 
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Figure 2. Radial velocity curves of V2083 Cyg for the primary 
(black points), and the secondary (red points), the lines represent 
the final fit (see the text and Tabic [2J. 



The visual orbit of the two components was derived by 
ISevmour et al.1 (J2002l ). They presented the orbital period of 
the double about 372 yr, the angular semimajor axis about 
498 mas, and the eccentricity 0.16. However, as they already 
mentioned, the orbit is still only a preliminary one. 



3 PHOTOMETRY AND SPECTROSCOPY 

We started collecting the photometric data of the system 
in April 2008. In total there are 31 nights of observations, 
but for the light-curve analysis we used only 27 nights of 
observations obtained from April 2008 to September 2009 
and carried out with the same telescope and detector at 
the private observatory by one of the authors (PS). Ow- 
ing to high brightness of the target, there was used only a 
small 34-mm refractor at the Private observatory in Brno, 
Czech Republic, using the SBIG ST-7XME CCD camera 
and s tandard BVR filters by the specification by iBesselll 
l|l990r) . All the measurements were processed by the soft- 
ware C-MuNlPACPO, which is based on ape rture photo metry 
and using the standard DaoPhot routines (|Todvll 19931 ). The 
other nights were used only for deriving the precise times of 
minima for the system. 

Besides our new observations, there were also used the 
photometric data ob tained within the SuperWASP survey 
IjPoUacco et al.ll2006r ). However, these data are not of enough 
quality to be used for the LC analysis. Hence, we made 
use of the SWASP photometry only for deriving the min- 
ima times of V2083 Cyg for a prospective period analy- 
sis. For al l of the minima the Kwee - van Woerden method 
was used (Kwee fc van Woerdenlll956l ). and all of them are 
given in Tabic lAll The linear ephemeris are as follows: 
HJD = 2448501.1237 + 1.867493429 ■ E. 

The CCD spectra were obtained at Ondi^ejov observa- 
tory, Czech Republic, using the 2.0-m telescope equipped 
with a SITe-005 800 x 2000 CCD detector. These spectra 
cover a wavelength region 626 - 676 nm. All of them were 
secured between April 2010 and May 2011 and have a lin- 
ear dispersion of about 17 A/mm. Their S/N values range 
typically between 100 and 200. 

For all of the spectra used the wavelength calibration 
was made via a ThAr comparison spectra obtained be- 
fore and after the stellar spectra itself. The flatfields were 
taken in the beginning and end of the night and their 
means were used for the reduction. After then, the radial 
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velocities (hereafter RV) were obtained with the program 
SPEFO (JHorn et al.lll996l ). with using the zero point cor- 
rection via measuring the telluric lines. In total 19 spectra 
were obtained this wa y. Moreover, also two Elodie spectra 
iJMoultaka et al.ll20o3 ) obtained in 1999 were added for the 
analysis. 

A list of derived radial velocities from all of the available 
spectral observations is written in Table [T] In the last col- 
umn the reference Elodie or Ondi'ejov is given. For all of the 
spectra we also tried to identify the third component lines, 
however these radial velocities are rather uncertain and af- 
fected by relatively large errors (see below Sections [S] and 

EJ. 



4 LC AND RV ANALYSIS 

The complete LC (in BVR filters) and RV curves were 
analyzed simultaneously, using the program PHOEBE 
(Prsa & Zwittcr 2005 ), which is based on th e Wilson- 
Devinney algorithm (|Wilson fc DevinnevI 119711 ). The de- 
rived quantities are as follows: the semi-major axis a, the 
mass ratio q = M2/M1, the systemic velocity 7, the sec- 
ondary temperature T2, the inclination i, the luminosities 
Li, the gravity darkening coefficients gi, the limb darken- 
ing coefficients Xi, the albedo coefficients Ai, and the syn- 
chronicity parameters Fi. The limb darkening was approxi- 
mated via linear cosine law, and the v alues of Xj were inte r- 
polated from van Hamme's tables, see Ivan Hammd (|l993). 

For the whole analysis, we followed this procedure: at 
the beginning, we fixed the temperature of the primary com- 
ponent at Ti = 7930 K (corresponding to spectral type A7, 
Cox 2000). We were trying to find the best LC-I-RV fit ac- 
cording to the lowest value of rms. There was a solution 
reached, but this one was unacceptable due to the fact that 
resulting values of Mi , M2 , Li , L2 , Ti , and T2 are in contra- 
diction with each other. In particular, the resulting spectral 
types as derived from AI, L and T differ significantly between 
each other. For this reason, we tried a different starting value 
of Ti . With this method we were changing the temperature 
Ti in the range from 8520 K to 7020 K (spectral types A3 to 
FO) and trying to find a consistent solution. For all of these 
attempts, the value of Ti remained fixed. 

Our final parameters as derived from the LC-I-RV fit 
are given in Table [51 The plot of the LC is shown in Fig. [TJ 
while the RV curves with the fits are given in Fig. [2] The 
value of eccentricity was fixed at 0. For the discussion about 
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Figure 3. Visual orbit of V2083 Cyg as displayed on the slty. The 
individual observations (black points) are connected with their 
theoretical positions (red points). The dashed line represents the 
line of the nodes, while the dotted one the line of the apsides. 
Parameters of the fit are given in Table [3l 



Table 1. Radial velocities of V2083 Cyg as derived from the 
spectra from the ELODIE archive, and from the Ondfejov 
observatory. 



HJD - 2400000 


RVl 


RV2 


RV3 


Ref. 


51405.3658 


95.594 


-81.861 




Elodie 


51407.4156 


129.872 


-117.06 




Elodie 


55316.5290 


20.920 


2.771 


-16.272 


OND 


55380.5210 


-122.394 


142.771 


-16.336 


OND 


55385.3570 


122.946 


-99.062 


-16.775 


OND 


55385.5530 


56.001 


-29.860 


-18.597 


OND 


55386.3640 


-73.970 


99.094 


-14.986 


OND 


55386.5650 


-.442 


25.184 


-17.640 


OND 


55405.3670 


55.349 


-26.402 


-16.327 


OND 


55425.3970 


-117.384 


141.598 


-16.145 


OND 


55496.2990 


-121.010 


140.992 


-16.273 


OND 


55496.4630 


-99.677 


124.485 


-16.277 


OND 


55497.2890 


141.451 


-119.863 


-15.212 


OND 


55622.5420 


116.769 


-96.782 


-15.048 


OND 


55622.6370 


89.474 


-66.160 


-14.284 


OND 


55671.4500 


-18.008 


36.024 


-14.574 


OND 


55671.5840 


-72.089 


92.471 


-14.950 


OND 


55671.6310 


-86.135 


105.183 


-14.784 


OND 


55689.5750 


144.545 


-122.291 


-16.380 


OND 


55692.3700 


-119.493 


143.646 


-15.231 


OND 


55692.5540 


-98.435 


122.422 


-15.080 


OND 



the physical parameters of the components (ecUpsing ones 
and also the third one), see Section [6l 

For the whole computation process, the values of albe- 
dos Ai and the values of gravity darkening coefficients 
were set at their appropriate values (Ai = 1 or 0.5, and 
Qi = 1 or 0.32) according to the component's temperatures 
(Ti < 7200 K, or T, > 7200 K). Another problematic issue 
were the values of Fi , which tended to decrease down to for 
both components for each of the Ti values. These dropped 
down very quickly after a few steps of iterations. Due to this 
reason we tried a different approach. From the spectra of the 
system we estimated the values vsini, which were used to 
derive the values of Fi for both components. Therefore, the 
values of Fi as given in Table [2] are not derived from the 
combined LC&RV analysis but from the spectra. 

The fitting process with PHOEBE was carried out as- 
suming three luminosities. Besides the luminosities of pri- 
mary and secondary component of the eclipsing binary pair 
also the additional third light La was considered. This lu- 
minosity corresponds to the visual component B and is pre- 
sented in the combined light for all the time (the two visual 
components are too close). From this value one can specu- 
late about some physical parameters of the third body in 
the system, see below Section [S] 



5 VISUAL ORBIT 

The close eclipsing pair is orbiting around a common 
barycenter with the third distant component of the system. 
Recent precise interferometric observations are to be used 
for determining the parameters of this visual orbit. Since 
its discovery as a double star by lAitkenI l|l904 ) there were 
obtained 61 astrometric observations of the double (i.e. po- 
sition angle and separation). We took these data from the 



Table 2. The final LC and RV parameters of V2083 Cyg. 



Parameter 


Value 


Parameter 


Value 


a[RQ] 


9.57 ± 0.15 


Li (B) [%] 


26.7 ± 0.7 


q = M2/M1 


0.97 ± 0.07 


L2 (B) [%] 


34.9 ± 0.9 


7 [km/s] 


10.78 ± 0.68 


L3 (B) [%] 


38.4 ± 0.8 


Ti [K] 


7630 (fixed) 


Li (V) [%] 


26.7 ± 0.6 


Ta [K] 


7623 ± 45 


L2 (V) [%] 


34.7 ± 0.9 


e 


(fixed) 


L3 (V) [%] 


38.6 ± 0.8 


i [deg] 


80.47 ± 1.60 


Li (R) [%] 


26.4 ± 0.6 


xi = X2 (B) 


0.412 


L2 (R) [%] 


34.2 ± 0.9 


xi = X2 (V) 


0.356 
0.356 


L3 (R) [%] 


39.4 ± 1.0 


XI = X2 (R) 


Derived physical quantities: 


91 =92 


1.000 (fixed) 


Rl [Rq] 


2.12 ± 0.17 


Ai=A2 


1.000 (fixed) 


R2 [Rq] 


2.45 ± 0.20 


Fi 


0.81 ± 0.13 


Ml [Mq] 


1.71 ± 0.11 


F2 


0.84 ± 0.11 


Ma [Mq] 


1.66 ± 0.09 



WDS database. The very last observation was obtained in 
2009. 

Since its discovery, the position angle of the pair 
changed of about 88 °. Thanks to this movement, the or- 
bit of the pair ar ound a barycen t er wa s derived. The orbit 
was published by ISeymour et al] l|2002l ). who computed its 
orbital period of about 372 yr. However, since this most re- 
cent study there were published three new interferometric 
observations, so we decided to perform a new analysis with 
the complete data set. 

Our new computation led to the visual orbit parameters 
given in Table [3] For the computation we used a following 
ap proach. Starting with the orbital parameters as published 
bv lSeymour et al.l (|2002l ). the final fit reached very different 
solution. Moreover, there were found several different min- 
ima in the parameter space as derived from this astrometric 
data set. Some minimum was found with very long orbital 
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Table 3. Final parameters of the long orbit. 



Parameter ^e^mour_et_aL (2002) This work 





^ 




" 


4.00 E, 


P3 [day] 


135869 


64778 ± 427 


4.20^ 

2 


P3 [yr] 


372 


177.4 ± 1.2 


4.40 £ 


a [mas] 


498 


291.9 ± 1.4 


4.60 


To 


2438395 


2400006 ± 375 


4.80 
5.00 


n [deg] 


73.6 


174.54 ± 2.9 


5.20 


a; [deg] 


189 


334.89 ± 5.3 




i [deg] 


64 


48.73 ± 3.6 




e 


0.16 


0.471 ± 0.018 



80000 100000 

Period p [d] 



Figure 4. Plot of period-mass and period-rms relation. The 
dashed lines represent our final fit. For details see the text. 



periods, but this solution seems to be less probable due to 
the poor coverage of the data. The most significant luinimum 
(the deepest one) was found near the period of 177.4 years. 
However, we would like to emphasize that the orbital so- 
lution is still a preliminary one. New precise observations 
secured every year would be very welcome for derivation of 
the orbital parameters with higher conclusiveness and es- 
pecially for setting the more solid constraints on pa and a 
values. These are the most important for a discussion about 
the nature of the third component (see Section [5] below). 

In Fig. |4] there is a plot of total mass versus pe- 
riod, as well as the rms of the particular fit versus pe- 
riod. For our final solution reached (minimum rms with 
P3 = 64778.357 day) the value of total mass was computed 
(using the Hipparcos parallax) - these are shown as dashed 
lines in Fig. 2] The relation between the two vertical axis 
(parallax and total mass) is defined via a third Kepler's law 
using our final solution. As one can see from the bold line of 
mass-period relation, the total mass as derived from our fi- 
nal solution is close to the minimal mass in this period range 
(the uncertainty of Hipparcos parallax tthip ~ 4.32 ±0.57 is 
shown as a grey area). Of course, this analysis is very sensi- 
tive to the input weightening scheme. The individual weights 
of the data points were set equal to each other, because for 
most of the observations the sigma or some other error es- 
timations are missing. N o minimum of rms ne ar a period of 
372 years as proposed bv lSevmour et al.l (|2002l ) is presented. 
One can ask, why such a different solution was reached when 
using only three new interferometric observations. The main 
reason (maybe besides different weightening) is that these 
three new measurements provide a strong constraints on the 
fit. This is due to the fact that the posit ion angle between 
our most recent data and those from Sc vmour et al.l (|2002l ) 
changed about 20°, which is about 1/4 of the total posi- 
tion angle range covered. All of these calculations (e.g. the 
Kepler's law) were using the set of reco mmended values of 
fundamental parameters as proposed bv lHarmanec fc Prsal 
l|201ll) . 

On the other hand, we also tried to compute the pre- 
dicted change of the third-body velocities over the time span 
of more than 11 years covered with our spectroscopic data. 
Taken into account some assumptions (masses), the change 



in velocity resulted in more than 20 km/s. Such a large 
velocity difference should be easily detectable in our RVs 
data. Unfortunately, we were not able to identify the third- 
component lines in the ELODIE spectra and in newer data 
from Ondfejov there is no such difference, hence we could 
only speculate about our findings. The reason could be ei- 
ther different masses or much longer orbital period. Another 
explanation is an incorrect identification of the third-body 
lines in the spectra. 



6 PHYSICAL PARAMETERS 

Taken into account all results as presented above, one can 
make a picture of the system, its geometry and orientation 
in space. From the combined LC and RV analysis there re- 
sulted that both eclipsing components are probably main 
sequence stars, located well within their respective Roche 
lobes. According to their masses and temperatures, it seems 
like their individual spectral types are probably of A7 and 
A8 (e.g. i_Poppor 1980, Harmanec 1988, or Andersen 199l|) 
for primary and secondary, respectively. However, according 
to their luminosities, it seems like the stars are of slightly 
earlier spectral type (about A5). 

Another task was to derive the value of the third light 
Ls from the LC solution and to obtain a magnitude differ- 
ence between the two visual components. This value resulted 
in about 0.49 mag, which is in rough agreement with the 
Am = 0.43 mag value as presented in the WDS catalogue. 

A discussion about the third body is still difficult due to 
some aspects of the problem. The most problematic issue is 
still the uncertainty of the Hipparcos value of parallax. The 
relatively high error of about 13 % could lead to distances in 
wide range from 204 to 267 pc. Thanks to this uncertainty 
also the value of total mass as computed from the visual 
orbit (see Table [3| lie in between 6.54 and 15.41 Mq with 
the mean value of 9.81 Mq. 

Subtracting masses of both eclipsing components, we 
obtain an interesting result of mass of the third body about 
6.44 Mq (with upper and lower limits about 12.24 and 2.97). 
Such a massive third body easily cannot be a main sequence 
A star as predicted from the Am value. One possible expla- 
nation of this discrepancy is that this component is also a 
double star. If we speculate about two identical stars, then 
such stars have to be of only slightly later spectral type than 
the eclipsing components, (because of the total mass). As- 
suming two FO stars, we can hardly satisfy the magnitude 
difference between the components. However, this explana- 
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tion is still questionable because the third lines in the spectra 
do not show a doubleness profile. 

To solve this discrepancy we also tried to use the pro- 
gram KOREL (Hadrava 2004 ) for disentangling of the spec- 
tra taken at Ondfejov observatory. However, it was also 
not able to solve the problem. The final parameters on 
one hand confirmed our findings about the LC-I-RV solu- 
tion (mass ratio q from KOREL resulted in about 0.993), 
but on the other hand it also results with the value of mass 
ratio qa = M3/M12 > 1. This would indicate that the third 
body is more massive than the eclipsing pair, but also the 
less luminous one. Solving the problem of its lower luminos- 
ity and higher mass with introducing a degenerate object is 
a highly speculative solution. Hence, the nature of the third 
body still remains an open question. The KOREL radial 
velocities of the third body were also used and these are the 
values presented in Table [1] in the RV3 column. 



7 DISCUSSION AND CONCLUSIONS 

The multiple system V2083 Cyg is still rather neglected one 
and this is the first detailed analysis of it. The components of 
an eclipsing binary are of spectral type A, are well-detached, 
with no evidence of circumstellar matter, emission in the 
spectra, etc. This close pair is also orbiting around a com- 
mon barycenter with the third component with period about 
177 years. The mutual inclination of the two orbits is 31.8°, 
therefore we can only speculate about a common origin of 
the system. 

A nature of the third component is still rather prob- 
lematic to derive. From the combined LC and RV analysis 
there results that the third body is slightly less luminous 
than the eclipsing pair. But the Hipparcos parallax indi- 
cates a higher total mass of the system than computed from 
all components' masses. A possible explanation is that the 
value of Hipparcos parallax is underestimated and a real 
distance of V2083 Cyg is lower (even outside of the error 
bars of the Hipparcos data). This would not be an excep- 
tional case, because for some systems the Hipparcos data 
yielded incorrect par allax due to the pre sence of a close vi- 
sual companion fe.g. iDocobo et al-lbOOST ). Another possible 
explanation is that the body is also a binary, but there are 
also some problems with this explanation (luminosity and 
the spectral lines of such body). For that reason, new more 
detailed observations are still very welcome. 

However, if the hypothesis of binarity of the third com- 
ponent is proved, it will shift the triple to quadruple. On one 
hand, such systems of higher multiplicity are of big interest, 
on the other hand we also deal with very incomplete statis- 
tics o f the m among the st ars (see e.g. lEggleton fc TokovininI 
l2008l . and lEggletonll2009l ). 
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Table Al. Heliocentric minima of V2083 Cyg 



HJD - 2400000 


Error 


Type 


Filter 


Observer 


54609.69354 


0.00179 


prim 


- 


SWASP 


54638.63794 


0.00035 


sec 


- 


SWASP 


54639.57227 


0.00042 


prim 


- 


SWASP 


54652.64554 


0.00089 


prim 


- 


SWASP 


54668.51720 


0.00046 


sec 


- 


SWASP 


54669.44747 


0.00046 


prim 


- 


SWASP 


54683.46156 


0.00027 


sec 


- 


SWASP 


54684.39467 


0.00095 


prim 


- 


SWASP 


54994.39887 


0.00091 


prim 


B 


PS 


54994.39888 


0.00113 


prim 


V 


PS 


54994.39948 


0.00128 


prim 


R 


PS 


55049.49292 


0.00110 


sec 


B 


PS 


55049.49015 


0.00104 


sec 


R 


PS 


55049.49055 


0.00062 


sec 


V 


PS 


55051.35585 


0.00081 


sec 


B 


PS 


55051.35567 


0.00078 


sec 


V 


PS 


55051.35645 


0.00107 


sec 


R 


PS 


55064.43157 


0.00091 


sec 


B 


PS 


55064.43162 


0.00057 


sec 


V 


PS 


55064.43102 


0.00064 


sec 


R 


PS 


55076.57570 


0.00116 


prim 


B 


PS 


55076.57197 


0.00108 


prim 


V 


PS 


55076.57194 


0.00165 


prim 


R 


PS 


55093.37749 


0.00087 


prim 


B 


PS 


55093.37848 


0.00068 


prim 


V 


PS 


55093.37624 


0.00063 


prim 


R 


PS 


55374.43431 


0.00021 


sec 


I 


RU 


55429.52450 


0.00028 


prim 


I 


RU 


55740.46289 


0.00056 


sec 


R 


PS 


55797.42284 


0.00043 


prim 


I 


RU 



Note: PS - Pctr Svoboda, RU - Robert Uhlaf 



